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Abstract 

The notion of a rough two-dimensional (convex) body is intro- 
duced, and to each rough body there is assigned a measure on T 3 
describing billiard scattering on the body. The main result is char- 
acterization of the set of measures generated by rough bodies. This 
result can be used to solve various problems of least aerodynamical 
resistance. 
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1 Definition of a rough set and statement 
of main theorem 

1.1 Introductory remarks and review of literature 

In this paper the notion of a rough two-dimensional (convex) body is given 
and some properties of rough bodies are established. 



"On leave from Department of Mathematics, Aveiro University, Aveiro 3810, Portugal 



1 



Let B C K 2 be a convex bounded set with nonempty interior, that is, 
a bounded convex body. Consider the "set" obtained from B by moving 
off a set of "very small" area. Such a (heuristically defined) set is called a 
rough body: from the "macroscopic" point of view, it almost coincides with 
B, and from the "microscopic" point of view, it contains some "flaws". (One 
can imagine a detail of a mechanism that, after a period of exploitation, has 
got some defects.) If the removed set adjoins the boundary dB, one can 
expect that a flow of point particles incident on the rough body is reflected 
in another way as compared to reflection from B. 

The notion of rough body arises naturally when studying Newton-like 
problems of the body of least resistance. The first problem of such kind 
was considered by Newton itself pQ. Recently there were made several works 
concerning the problem of least resistance in various classes of admissible 
bodies; see, e.g., [2]-[T3], [15]. The solution of a minimization problem for 
the case of rotating bodies can be naturally identified with a rough body ([H]; 
see also concluding remarks to this paper). 

There are many papers on particle scattering by rough bodies (see, e.g., 
[IS]- [IB]); they describe bodies and flows of particles that occur in nature. 
On the contrary, we assume that a rough structure can be "manufactured", 
and our aim is to describe all possible rough structures. 

1.2 Definition of a rough body 

It is supposed that the "microscopic structure" of the boundary of a rough 
body can be detected from observations of particle scattering on the body. 
From this point of view, two rough bodies are considered equal if they scatter 
flows of particles in an identical manner. Having these observations in mind, 
we give the definition of a rough body. 

Let B be a bounded convex body. Denote by n(£) the unit outer normal 
vector to dB at a regular point £ G dB, and denote by (dB x S 1 )-^ the 
set of pairs (£, v) G dB x S 1 such that (n(£), v) > 0. Here and in what 
follows, ( •, •) means the standard scalar product in IR 2 . The set (dB x S l ) + 
is equipped with the measure \x which is defined by v) = (n(£), v) d£ dv, 
where d£ and dv are the one- dimensional Lebesgue measures on dB and S , 
respectively. 

Let Q be a set with piecewise smooth boundary contained in B\ consider 
the billiard in IR 2 \ Q. Note that Q is not necessarily connected. For (£, v) G 
(dB x S 1 ) + , consider a billiard particle starting at the point £ with the 
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velocity —v. After several (maybe none) reflections from dQ\dB, the particle 
will intersect dB again, at a point £ + = £q B (£,v) G dB; denote by v + = 
Vq B (^,v) the velocity at this point. It may happen that the initial point £ 
belongs to dQ; in that case we have £ + = £ and the vector v + is symmetric to 
v with respect to It may also happen that at some moment the particle 
either gets into a singular point of dQ, or touches dQ at a regular point, or 
stays in B \ Q forever and does not intersect dB again, or makes an infinite 
number of reflections in finite time. The set of corresponding points (£,i>) 
has zero measure, and the corresponding values £q B (£,v) and-y^ B (£, t>) are 
not defined. 




Thus, there is defined the one-to-one mapping Tq^b '■ (£, v) \— > (£q i>), f q f )) 
of a full measure subset of (dB x S' 1 ) + onto itself. It has the following prop- 
erties: 

Tl. Tq B preserves the measure \i. 
T2 - T q!b = T Q,b- 

The mapping Tq B induces the measure vq B on T 3 = 5 1 ' x 5 1 x S 1 in the 
following way. Let A C T 3 be a Borel set; by definition, 

^ QiB (A) = /i ({(£, V ) G (dB x ^) + : ( W , v+ <B (Z,v),n(0) G A}) . 

In fact, the measure z/q^ contains information about particle scattering on 
Q. Imagine that an observer has no means to track the trajectory of particles 
inside B. Instead, for each incident particle there is registered the triple of 
vectors: the initial and final velocities (measured at the points of first and 
second intersection with dB), and the normal vector to dB at the point of first 
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intersection with dB. The normal vector at the second point of intersection 
is not registered; as will be seen later on (lemma 1), if the area of B \ Q is 
small then the difference between the normal vectors at these two points is 
also small. The measure vq,b describes the distribution of triples. 

Definition 1. We say that a sequence of sets {Q m , m — 1, 2, . . .} represents 
a rough body, if 

Ml. Q m C B and Area{B \ Q m ) ^ as m — > oo; 
M2. the sequence of measures VQ m ,B weakly converges. 

Two sequences of such sets are called equivalent, if the corresponding limiting 
measures coincide. An equivalence class is called a body obtained by rough- 
ening B , or simply rough body, and denoted by B, and the corresponding 
limiting measure is denoted by u B . 

Note that the sets Q m in this definition are not necessarily connected. 

Remark. Since T 3 is compact and the full measure o/T 3 satisfies z/q^TT 3 ) < 
27T \dB\, one concludes that the set of measures {vq,b}, with fixed B , is weakly 
precompact. That is, any sequence of measures {vQ mt B} contains a weakly 
converging subsequence. In this sense one can say that a sequence, satisfying 
only the condition Ml, can represent more than one rough body. 

We would also like to mention that, firstly, two rough bodies obtained 
one from another by translation are identified, according to our definition. 
Secondly, particle scattering on E in a small neighborhood of £ G OB can 
be detected if £ is an extreme point of B, and cannot otherwise. Indeed, 
if £ is an extreme point of B, the scattering is described by the restriction 
of vb on T 2 x J\f n (£), with A/" n (£) being a small neighborhood of n(£) in S 1 . 
If, otherwise, £ is not an extreme point of B, that is, belongs to an open 
linear segment contained in dB, the scattering can only be determined on 
the whole segment. 

Actually, from the viewpoint of applications to the problems of optimal 
resistance in homogeneous and rarefied media (see the section Concluding 
Remarks and Applications), these drawbacks are not so serious. Indeed, 
resistance of a body is invariant under translations (due to homogeneity). 
Besides, if the boundary of a body contains a linear segment, one does not 
need to know scattering at each point of the segment; it suffices to know it 
on the segment in the whole (due to homogeneity and rarefaction). 
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The definition of a rough body could be made in a slightly different way 
basing on measures defined on 5 1 x 5 1 x dB. In that case the triple (v, f + , £) 
should be registered, with £ being the point of first intersection with dB. 
That definition would allow one to register particle scattering at each point 
of dB and to distinguish between bodies obtained by translation one from 
another. However, we prefer to adopt the former definition, since it seems to 
us mathematically more transparent and makes the arguments a bit easier. 

1.3 Examples 

Sometimes it is convenient to use another representation of the measure v&. 
Namely, consider the change of coordinates (v,v + ,n) i— > (tp, tp + , n), where 
tp = Argf — Argn, <p + = Argt> + — Argn. Here Argt> is the angle be- 
tween a fixed vector and v measured, say, clockwise from this vector to v. 
If (v,v + ,n) G spti^B then (p and tp + belong to [— 7r/2, it/2] modulo 2ir. In- 
troduce the shorthand notation □ := [— n/2, 7r/2] x [— 7r/2, 7r/2] and de- 
fine the mapping w : □ x S* 1 — > T 3 by w(tp,tp + ,n) = (v,v + ,n). One has 
spt ujb C zu(n x S* 1 ). Denote pq := (ro^ 1 )*//^. Sometimes this measure can 
be factorized: uq = tjb <S> Tb, where r/s is defined on □ and tb is the surface 
measure on B; so to say, the "roughness" is "homogeneous" along the body's 
boundary. Consider several examples. 

Example 1: "smooth body". 
The rough body represented by the sequence Q m = B is identified with B 
itself. The corresponding measure is vb — Vo® t b, where the measure 770 h as 
the density cosy? • 5((p + ip + ); the support of i] is shown on the figure below. 
On this figure, B is taken to be an ellipse. 




,+ 
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Example 2: roughness formed by triangular hollows. 
Q rn is a 2m-polygon; the angles 270° alternate with the angles that are 
slightly smaller than 90°. All vertices corresponding to the angles smaller 
than 90° belong to dB. Any two sides that form an angle 270° are equal. 
The largest side length tends to zero as m — > oo. Thus, the set Q m is obtained 
by moving off m "hollows" from its convex hull, each of the hollows being an 
isosceles right triangle. 

The corresponding measure is vq = r] v <g> tb, where the measure ^ v has 
the density 



cos tp 



7I\ 



X[-*/2,-*/4\((p) 6(<p + <P + + ^) + X[-*/4,*/*](<P) $(<P - <P + ) + 



TT , 



+X[n/4,n/2](<p)S((fi + if + - - 



+ | sin 



7I\ 



X[-*/4,0](<P) 8(<P + <P + o)- 



-X[-n/A,n/4]((p)S((p - ¥ + ) + X[0,7r/4](^) 5{ip + ^ ~ 



Thus, the support of i] v is the union of three segments; see the figure below. 
The middle segment y? + = tp corresponds to double reflections, and the lat- 
eral segments, y? + = — ip — tt/2 and <p + = —tp + tt/2, correspond to single 
reflections, from the right or from the left side of a triangular hollow. On the 
figure, B is a circle. 





Example 3: roughness formed by rectangular hollows. 
The sets Q m are obtained by removing a finite number of "rectangular hol- 
lows" from B. In other words, one has Q m = B\(U n Q mjn ) , where the removed 
sets Q m , n do not mutually intersect and each set dQ m ,n \ dB is the union of 
three sides of a rectangle. The ratio (width)/(depth) of a hollow depends 
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only on m and is denoted by h m . Denote by l m = \dB \ U n (dfl m<n ) \/\dB\ 
the relative length of the part of boundary OB not covered by hollows. We 
assume that linim^oo h m = = lim m ^oo l m . On the figure below, B is a 
square. 



ruuuuuuuuuwuuuun 
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The measure equals v B = <S> t b . The density of the measure r/n 
equals \ cosy? ■ (5((p + tp + ) + 5(ip — y? + )), and the support is the union of two 
diagonals, (p + = if and ip + = —ip; see the figure. The particles with even 
(odd) number of reflections contribute to the first (second) diagonal. 



1.4 Main theorem 

According to the definition (TJ each rough body is identified with a measure 
on T 3 . The question is: what is the set of these measures? The following 
definition and theorem give the answer. 

Let us first introduce some notation: n v ^ n : T 3 — > T 2 , n n : T 3 — >• S 1 , 
etc. are projections onto the corresponding subspaces: TT Vin (v, v + , n) = (v,n), 
7r n (t>, v + , n) = n, etc.; n d : T 3 — > T 3 is the symmetry with respect to the plane 
v = v + , that is, 7Td(v, v + , n) = (v + ,v,n); z + = max{0,^} is the positive part 
of z G R; and u means Lebesgue measure on S 1 . Recall that tb is the surface 
measure on B and is defined on S 1 . 

Definition 2. We denote by Mb the set of measures v on T 3 such that 
Al the marginal measures vr* n z/ and nf + u are 

A2 7rf v = v. 



7 



Denote also M. = Ub-M-b, the union being taken over all bounded convex 
bodies B. 

Taking into account the Alexandrov theorem on characterization of sur- 
face measures, one concludes that M. is the set of measures v on T 3 such 
that 

1) the marginal measure 7r*z/ =: r satisfies the conditions 

la. f sl n dr{n) = 0; 

lb. for any v G S 1 holds f sl (n,v) 2 dr(n) ^ 0; 

2) the marginal measures n# n v and nf + n v satisfy the conditions 

2a. TT# n u = (v,n) + -u®T; 

2b. nf+^v = (v + , n)+ -u®t. 
Thus, these marginal measures coincide; the only difference is in the no- 
tation for the variables: v, n in the case 2a and v + ,n in the case 2b. 

Now we can state the main theorem. 

Theorem. The set of measures {vb}, with B being all possible bodies 
obtained by roughening B, coincides with M.b- Therefore, {z/g, B is a rough 
body} = M. 

In section 2, we formulate two auxiliary lemmas and using them, prove the 
theorem. In section 3, the lemmas are proved. Section 4 contains concluding 
remarks and applications of theorem to problems of optimal aerodynamic 
resistance. Appendices A and B contain proofs of some auxiliary technical 
results. 



2 Statement of auxiliary lemmas 
and proof of theorem 

2.1 Statement of lemma 1 

Fix a bounded convex body B. Two points £i, £ 2 G dB, £ x ^ £ 2 divide the 
curve dB into two arcs. Denote by £2) the length of the smallest arc and 
denote 

. f 16-61 

c = c B := mf — ; 
one obviously has < c < 1. 
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Let Q C B; denote 



l£-£ + l Q , B := // 

<(dBxS 1 )+ 



and 



n — rv 



H0-n(& B (£,v))\d^,v). 
Lemma 1. (a) The following holds true: 



IB — / / rns; '"\^q,i 

'(9Bx5 1 )^ 



|e-e + l g , fl <27r-Area(S\Q). 
(b) For sufficiently small Area(£> \ Q) .0 one has 

27r\/87r 



n — n~ 



\ Q)B < — VArea(5\g). 



2.2 Statement of lemma 2 

Let us first introduce the notion of a hollow. 

Definition 3. Let Q G M 2 be a closed bounded set with piecewise smooth 
boundary and I C dQ, where 

(i) I is an interval contained in a straight line (x,n) = a and 

(ii) Q\I is contained in the open half-plane (x,n) < a. Here n is a fixed 
unit vector. 

Then the pair (Q, I) is called a hollow oriented by n, or just an n-hollow. 

Here and in what follows, I is shown dashed, and dVt \ I is shown by solid 
line. 



That is, smaller that a value depending only on B. 
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(x, n) = a 



Define the measure jli on I x S* 1 by djli(^,v) = m + d£ dv, where 
|/| means the length of /. Obviously, flj is supported on the set (/ x 
S 1 ) + := {(£,v) E I x S 1 : (n,v) > 0}. Define the one-to-one mapping 
(£, v ) l— f), ^"j^, u)) of a full measure subset of (I x S* 1 ^ onto 
itself. Namely, consider the billiard in Q. Let (£, i>) G (/ x S' 1 ) + ; consider 
the billiard particle starting at the point £ with the velocity —v. It makes 
several reflections from dfl \ I and then reflects from / again, at a point 
S + = E^j(£,v). The velocity immediately before this reflection is denoted 
by V + = Vq j(£, v). The mapping so defined preserves the measure jli and is 
an involution, that is, coincides with its inverse. 

One can give an equivalent definition based on the mapping C,q B (C,, v), vX B (C,, v) 
just defined in subsection 1.2. Take a set Q such that Q is a connected 
component of convQ \ Q and / is a connected component of <9(convQ) \ 
dQ. For (£,v) E (I x S 1 )*, let by definition (E+^,v), V^fov)) : = 

(^Q,convQ(^ w )> v Q,convQ(^ v ))- This definition does not depend on the choice 
of Q. 

Definition 4. Let (Q, 7) fre a hollow. The measure rjnj on T 2 = S 1 x S* 1 
defined as follows. For a Borel set AcT 2 , put 

rjn,M) :=i*i({(M S') + : (v, V^fot;)) G A}). 

We shall say that rjnj is the measure generated by the hollow (£1,1). 

Here we use the notation ir v , ir v + : T 2 — > S* 1 for the projections onto 
the subspaces {v} and {t> + }, respectively; n v (v,v + ) = v, n v +(v,v + ) = v + . 
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We also denote by iid the symmetry with respect to the diagonal v = v + ; 
n d (v,v + ) = (v + ,v). 

Definition 5. Denote by A n the set of measures r\ on T 2 such that 

1) d,7rfr](v) = (v,n) + dv, dixf + ri{v + ) = (v + ,n) + dv + ; 

2) Tx*n = 77. 

Any measure rjnj generated by an n-hollow belongs to A n . Indeed, for 
any A C S 1 one has n#r] n j(A) = mj {A x S l ) = v) E (/ x S r ) + : v E 

A}) = pi JJ IxA (n,v) + dC, dv = J A {n,v) + dv. This proves the first equality in 
1). 

Similarly, one has n*+r]n,i{A) = rjnjiS 1 x A) = E (I x : 

Vqj(C,,v) E A}). Since the mapping (£,i>) 1— > (Sq 7 , V^~ 7 ) preserves the mea- 
sure, one gets the value E (I x S' 1 ) + : t> G A}), which in turns 
equals to J A {n, v) + dv. This proves the second equality in 1). Finally, the re- 
lation 2) for rjaj is a simple consequence of involutive and measure preserving 
properties of the mapping (£,i>) 1— > (S^ 7 , Vq~ 7 ). 

Lemma 2. TTie set 0/ measures generated by n-hollows is weakly dense in 



2.3 Proof of the direct statement of theorem 

Here we prove that for any body B obtained by roughening B holds v& E Mb- 
Let Q C B; define the measure Vq B on T 3 by 

v^ B {A) := p ({(£,„) E (dB x S r ) + : („, G A}) , 

where A is an arbitrary Borel subset of T 3 . Thus, the definition of both vq B 
and Uq B is based on observations of vector triples (v,v + ,n) and (v,v + ,n + ), 
respectively. Here n and n + are the outer normals to dB at the points where 
the particle gets in B and gets out of B. The measures vq B and Vq B have 
the following properties: 

^t,n v Q,B = (v,n}+-u®T B , (1) 
*t,n + VQ,B = (v + ,n + ) + -u®r B , (2) 
^t v Q,B = Vq jB - ( 3 ) 
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Consider a sequence {Q m } representing £>; let us show that ^Q m ,B~^Q mi B 
weakly converges to zero as m -> oo. It is enough to prove that for any 
continuous function / on T 3 holds 

f(v, v + , n) dv QmiB {v, v + , n) - / f(v, v + , n + ) dv' Q B (v, v + , n + ) -^ m _»oo 0. 

(4) 

Taking into account the formulas for change of variables 



f(v, v + , n) du Q!B (v, v + , n) = / f(v, vl B (£, n(0) , v) 

T3 J(dBxS 1 )+ 



and 



f(v, v + , n + ) dv' QB (v, v + , n + ) = / f(v, v% B (£, v), 7i(£o,fl(f > v))) dfi(£, v), 

T 3 ' J(dBxS 1 )+ 

the formula (ED takes the form 



" wo ° J(dBxS 1 )+ 

(5) 

According to lemma (H the difference w(£g mB (£, v)) — n(£) converges to 
zero in mean, hence it converges to zero in measure; therefore the difference 

also converges to zero in measure. It follows that the formula §5§ is true. 

Thus, both VQ m B and v'n mB weakly converge to v&. Substituting Q = Q m 
into the formulas ([IH3]) and passing to limit as m — > oo, one gets 

^tn^B = (v,n) + -U(g)T B , 

7T* = (v + , n)+ • u ® re, 



that is, z/g G .Mb- 
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2.4 Proof of the inverse statement of theorem 



Here it is proved that for any v G M. b there exists a body B obtained by 
roughening B such that vjs = v. The proof is based on two statements. 

Statement 1. Let B be a convex polygon. Then for any measure v G A4b 
there exists a body B obtained by roughening B such that z/g = v. 

Proof. Let us enumerate the sides of the polygon B and denote by the 
length of the ith side, and by n i: the outer unit normal to this side. By 8 n , 
denote the probabilistic atomic measure on S 1 concentrated at n G S 1 , that 
is, 5 n {n) = 1. The surface measure of B is r B = ^Cj<5 ni ; this implies that 
any measure v G Mb has the form v = Yl c iVi ® S nv where rji G A nv 

According to lemma [2X any measure rji is the weak limit as m — > oo of 
measures rjnrnj™ generated by a sequence of rij-hollows (fi™, If 1 ). Now take a 
sequence of sets Q m such that conv Q m = B and each connected component 
of B \ Q m is the image of a set f2 4 m under the composition of a homothety 
with positive ratio and a translation, and additionally, the image of i] 71 under 
this transformation belongs to (ith side of B) \ dQ m . We also require that 
Area(£> \ Q m ) — > and | (ith. side of B) \ dQ m \ =: c™ — > Q as m — > oo. 
In Appendix A it is shown how to construct such a sequence Q m . The 
measure VQ m B = v m + v™ is the sum of the measure v m corresponding to 
reflections from OB n dQ m and the measures u™ corresponding to particles 
getting into the "hollows on the ith side" . One has u m = Yi( c i~ c T) -r ?o® ^ 
and ^ m = c™ • 7)si™,i™ ® ^ n i- ^ ne norm °f goes to zero and v™ weakly 
converges to Ci rji ® <5 ni for any i; it follows that VQ m ,B weakly converges to 
z/ as m — > oo. Therefore, the sequence Q m represents a body B obtained by 
roughening B, and Pq = v. □ 

Statement 2. For any measure v G A4b t/iere exzsi a sequence of convex 
polygons Bk C B with Area(B \ Bk) — > and a sequence of measures z/& G 
A^B fc weakly converging to v as k — ► oo. 

Proof. Consider a partition of the circumference S 11 into a finite number of 
arcs, S 1 = Uj«S l . It induces the partition of dB into arcs 95* = {£ G (95 : 
G iS 1 }. Consider the polygon B inscribed into dB whose vertices are 
separation points of this partition. Denote by the outer normal to the ith 
side of this polygon. Denote by s VljV2 the operator of rotation on S 1 that 
takes Vi to t>2, and define the mapping Tj : T 2 x S" 1 — > T 2 by Ti(v,v + ,n) = 
( s n,m v , s n, ni v + )- Finally, consider the measure v = YliWlVi ® $nv where 
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I b l \ is the length of the ith side of the polygon, and the measure rji on T 2 
is defined by T]i(A) = r^nv(T^ (A)) for arbitrary Borel set A C T 2 . Here 

\dB l \ is the length of the arc OB' 1 . One easily verifies that v belongs to M. B . 

Now take a sequence of partitions of S 1 , {<S k }i, k = 1, 2, . . ., where the 
maximum arc length of a partition goes to zero as k — > oo. Denote by 
{<9£>^}i, = 1, 2, . . . the sequence of induced partitions of dB, and take 
the sequence of polygons B k generated by these partitions. One clearly has 
Area(B \ B k ) -> and 

16* I 

max , - k . , — > 1 as A; — > oo, (6) 
i \dB\\ V ; 

where \b\\ is the length of the ith side of B k . In the same way as above, one 
defines the mappings T ik : T 2 x S l k — ► T 2 and the measures ^ = |6J.| 77^ eg) 
5 njfc G M Bk , where is given by r) ik (A) := z/(T^(A)) and n ik is the 
outer unit normal to the ith side of B k . 

It remains to show that v k weakly converges to v. For any continuous 
function / on T 3 one has 



f(v } v + ,n)du k (v,v + ,n) = ^2\b l k \ // f(v,v + ,n ik )dri ik (v,v + ) = 

f(T ik (v, v + , n),n ik ) du(v, v + , n). (7) 



T 3 ; JJ^ 2 

b 



E 



dB 



For each k define the mapping from T 3 to T 3 by the relations (v,v + ,n) 1— > 
(Tik(v, v + , n), n ik ) if n G S\. It uniformly converges to the identity mapping 
as k — » 00; hence the function /fc, defined by the relations f k (v,v + ,n) : = 
f(Y ik (v, v + , n), n ik ) if n G <S|. , uniformly converges to / as — * 00. From 
here and from ([6]) it follows that the right hand side in (J7D converges to 
iU T 3 f( v , v+ , n ) dv(v, v + ,n) as — > 00. Thus, the convergence J jdv k — > 
J fdu is proved. Q.E.D. □ 

The inverse statement of the theorem follows from statements 1 and 2. 
Indeed, let v G M.b- Using statement [2l find a sequence of convex polygons 
B k C B and a sequence v k G -M_B fc weakly converging to v. According to 
statement [fl each measure v k is generated by a rough body. Consider the 
sequence of sets Q k i C B k , I = 1, 2, . . . representing this body, and then 
from all these sequences choose a diagonal sequence Q k = Qu k such that 
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the corresponding sequence of measures VQ k B weakly converges to v and 
Area(.B \ Qk) goes to zero as k — > oo. The sequence Qk represents a body B 
obtained by roughening B and vs — v - 

3 Proof of the lemmas 
3.1 Proof of lemma [I] 

Consider the billiard in IR 2 \ Q. For (£,v) G (dB x S* 1 )-^, denote by r(£, f) 
the time the billiard trajectory with the initial data £, — v spends in B \ Q. 
In particular, if £ G dB fl dQ, one has r(£, v) = 0. 

Denote by D the set of points (x, w) G (B\Q) x S 1 that are accessible from 
(dB x S 1 )-^; that is, there exists G (95 x S' 1 ) + such that the billiard 

particle with the data £, — f at the zero moment of time, at some moment < 
t < t(£,v) will pass through x with the velocity w. This description defines 
the change of coordinates in D : (£,v,t) h- > (x,w); (£,f) G (dB x S* 1 )^ t G 
[0, t(£,v)], and the element of phase volume d 2 xdw in the new coordinates 
takes the form dfi(£, v) dt. Hence, the phase volume of D equals JJf D d 2 x dw = 
ffr gBxS i\ T~(£,v)d/j l (£,v). Taking into account that D C (B \ Q) x S 1 and 
the phase volume of (B\Q) x S 1 equals 2tt ■ Aiea(B \ Q), one gets 



This is in fact a simple modification of the well-known mean free path formula 



One has r(£, v) > |£ — £q B (£, v)\: the time the particle spends in B \ Q 
exceeds the distance between the initial and final points of the trajectory. 
This inequality and (jSJ) imply (a). 

The points £ and £,n B (£,,v) divide the curve dB into two arcs; denote 
by l(£,i v ) the shortest one. One has |7(£, i>)| = l(£, Co v ))i therefore 
\C-Cq,b(C,v)\>c\j(C,v)\. It follows that 




(8) 



(see, e.g., [TJ]). 




(OBxS 1 ) 



(dBxS 1 ) 



(9) 
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Let g(y) be a natural parametrization of the curve dB, g : [0, \dB\] — » 
95. By /(y) denote the measure of the values (£, i>) such that the in- 
terval 7(£, i>) contains the point g(y); that is, /(y) := jj(Q BxS 1 ) + ^(@(y') e 
7(£,t>)) dfj,(^,v). Making change of variables in the integral in the left hand 
side of (jOJ), one gets 



// 

J J (i 



r\9B\ 

\n/(€,v)\d(i(£,v) = / f(y)dy, 

(8BxS 1 ) + JO 



therefore 



r\9B\ 2 
/ f{y)dy<—Axe&(B\Q). (10) 

One easily sees that |/(y x ) - /(y 2 )| < 4 |yi - y 2 | for any yi and y 2 
and /(y) > 0. From here and from (ITU]) it follows that for sufficiently 
s mall A r ea(5 \ Q) (na mely for Area(5 \ Q) < c\dB\ 2 /(2ir)) holds f(y) < 
V / 87r/cv / Area( J B\Q). 

Recall that Arg(-u) is the angle the vector v ^ forms with a fixed vector 
vq] the angle is measured clockwise from vq to v and is defined modulo 27r. 
Introduce the shorthand notation £ + := £q B (£,v) and denote by AArg(£,t>) 
the smallest in modulus of the values Arg(n(£ + ))— Arg(n(£)). In other words, 
AArg(£, v) equals to the smallest of the values 

/ |dArg(7i c 0l, / |dArg(n € OI- 

J^,v) JdB\y(£,v) 

Taking into account that |n(£ + ) — n(£)| < |AArg(£, v)\, one gets that 



K£ + )-n(£)|< / |dArg(n^ 



and therefore, 



\ n ~ n+ \ Q ,B < // (/ \dArg(n^)\) dfi(^v) 

JJ(dBxs i )+ \V 7 (e,«) / 

Making change of variables in this integral, one obtains 

-\8B\ 



\n-n+\ QB < 



\ f(y) |dArg(n e(tf) )| < 27T V ^7^ v /Area( J B\Q). 
Jo 



Thus, (b) is also proved. 
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3.2 Proof of lemma [2] 

Fix n G S 1 and m G N. Let a be an involutive permutation of {1, . . . , m}, 
that is, cr 2 = id. Divide the half-circumference := {v G S* 1 : (t> , n) > 0} 
into m arcs S^ m = S„, . . . , 5™ m = 5™ numbered clockwise, such that for 
any i, f qi (v,n) dv = 2/m. For the sake of brevity we omit the subscript m 
when no confusion can arise. 

Definition 6. A measure n is called a (a, n) -measure i/i)6A„ and sptrj C 
U™ 1 x Sn^ % j, and therefore, for any i holds rj (^S^ x Sn^^j = 2/m. 

Proposition 1. For any measure rj G A n there exists a sequence of involutive 
permutations on {1, ... , m^}, k = 1, 2, . . . snc/i t/iat m^ tends to infinity 
and any sequence of (o ~k,n) -measures weakly converges to n as k — > oo. 

Proposition 2. Lei cr 6e an involutive permutation on {1, . . . , m}. Then the 
distance (in variation) between the set of measures generated by n-hollows and 
the set of (a, n) -measures does not exceed 16/m. In other words, whatever 
e > 0, there exist a (cr, n)-measure rj and an n-hollow (Q, I) such that \\t]qj — 
rj\\ < 16/m + e; here the norm means variation of measure. 

This distance actually equals zero, but we only need the (weaker) claim 
of proposition El 

Lemma [2] follows from propositions [1] and [2j Indeed, let rj G A„. First 
choose the sequence of permutations cr^, according to proposition [T], and then, 
using proposition [21 for every k choose an an n-hollow (ilk, Ij.) such that the 
distance from r)n hl i h to the set of (o^, n)-measures does not exceed 17 /m k . 
The sequence of chosen measures i]n k ,i k weakly converges to n. 

3.3 Proof of proposition [I] 

Introduce on the angular coordinate <p = Argv — Argn; that is, <p changes 
between — 7r/2 and 7r/2 and increases clockwise. With this notation, to the 
arcs S l n m correspond the segments J l m = [arcsin(— l+2(z — l)/m), arcsin(— 1 + 
2i/m)]. Define the measure A on [— tt/2, tt/2] by d\(<p) = cos (p dtp and de- 
note by A the set of measures n on □ := [—tt/2, tt/2] x [—tt/2, tt/2] such 
that (a) 7r*n = A = ir^+rj and (b) 7r*n = n. Here 7r v , 7r v +, and ir d 
are defined by ir^ip,^) = ip, 7r lp +(ip,ip + ) = ip + , ir d (ip,<p + ) = {<p + ,p). Re- 
formulating definition El we shall say that n is a a-measure if n G A and 
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spt ?7 C U™ x (j^ x J™ (i) ). Notice that in the new notation, the objects do 

not depend on n anymore: we write A instead of A n , a- measure instead of 
(a, n)-measure, and hollow in place of n-hollow. 

In this notation, proposition [T] can be reformulated as follows: for any 
measure rj e A there exists a sequence of involutive permutations a k on 
{1, . . . , m k }, k = 1, 2, . . . such that mk tends to infinity and any sequence of 
CTfc-measures weakly converges to r\ as k — > oo. 

The idea of the proof is as follows. First, r] is approximated by means 
of a rational matrix, and then, this matrix is approximated by means of a 
larger matrix generated by a permutation. 

Consider the partition of □ into smaller rectangles D k J = J l k x J k , i, j — 
1, . . . , k. Choose rational nonnegative numbers c k such that c k = c J k , ^2 ■ c k = 
2/k for any i, and \r/ (D^) — c l k 3 \ < k~ 3 for any i and j. To do so, it suffices 
to take positive rational values c k such that r](O k ) — k~ A < c k < r]{D k ) 
for i > j and put = c£ for i < j and = 2/k — J2j^i c k f° r ^ = 3- 
One has r? (Jjj x [-vr/2, tt/2]) = £j =1 77(D£) = 2/ A;, hence cf - r/(Df) = 
Ej/< («) - cjf) e [0, (* - 1) • A;" 4 ] C [o, k-% 

Any sequence of measures r\ k satisfying the conditions r] k {D k ) = c k , 
1 < j < A: weakly converges to rj. Indeed, for any continuous function / 
on □ holds 

Jo Jo fj^x H 

as k — > oo. 

To complete the proof, it suffices to find an integer m k > k and an invo- 
lutive permutation a k of {1, . . . ,m k } such that any cr^-measure, r] k , satisfies 
the equalities r] k (D k ) = c k , i, j = l,...,k. Choose a positive integer N 
such that all the values a^- := N ■ c k are integer. The obtained matrix 
A = (ciij)ij =l is symmetric and for any i the value Y^j=i a ij = ^N/k is a 
fixed positive integer. In Appendix B it is shown that there exist square ma- 
trices Bij = (bij)^v of size 2N/k such that Bj- = Bj i} the sum of elements 
in any matrix B^ equals and the block matrix D = (B^) composed of 
these matrices has exactly one unit in each row and each column, and other 
elements are zeros. 

D is a symmetric square matrix of size 2N; denote its elements by dy. 
Put m k = 2N and define the mapping o k on {!,..., 2A^} in such a way that 



18 



dia k (i) = 1 for any i. The so defined mapping o k is a permutation; it is 
involutive since the matrix D is symmetric. Moreover, if % is a o^-measure 
then for any i and j holds ^(D^ ) = A^" 1 ^ v = c£ . The proposition is 
proved. 

3.4 Proof of proposition [2] 

1. Whatever the n-hollow (Q, I), one introduces the reference system (xi, £2) 
in such a way that n coincides with (0, —1), and the interval / belongs to 
the straight line x<i = and contains the origin O = (0, 0). Like in the proof 
of proposition [Tj, introduce the coordinate p = Argt> — Argn on S^. One 
has v = — (simp, cosy?), <p G [— 7r/2, 7r/2]. The definition of the segments 
J % m = J\ the measure A, the set of measures A, and the cx-measure see in 
the beginning of the previous subsection. The mapping (£, v ) 1— > Vqj(£, v) 
in the new coordinates £, ip is written as (£, <p) t—> p^j(^,p). Finally, define 
the measure /i/ on / x [— n/2, ir/2] by <i/i/(£, </?) = <i£ <iy9. 

Denote □' = (u™^) x (u™^) , = J 1 x [-tt/2, tt/2], D 2 = 
J m x [-tt/2, tt/2], D 3 = (UJS 1 J*) x J 1 , and D 4 = (U^J*) x J m . Thus, 
one has □ \ □' = Dx U D 2 U D 3 U D 4 ; see the figure below. 





□4 




□1 


□' 


□2 




□3 





It suffices to construct a sequence of hollows (fi e , J e ), e > such that 

(P) for any i / 1, m, cr(l), cr(m) i/ie measure of the set of 
values (£,,p) G I E x J* swc/i £/ia£ ^ 7 ((, ^ goes to zero 
as e — ► 0. 

Then, speaking of restrictions of measures on the subset □', one gets that 
the distance from the restrictions of measures rja e j e to the set of restrictions 
of a- measures goes to zero as e — > 0. On the other hand, for any measure 
77 G A one has ^(Di) = ^(Cta) = 2/m, ^(Ds) < 2/m, 77(^4) < 2/m, hence 
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r](n \ □') < 8/m, therefore the distance between the restrictions on □ \ □' 
of any two measures r/i and r] 2 from A does not exceed 16/m: || ^iJ D \ D ' — 
V2\n\D' II — 16/m. It follows that the upper limit of distances from rjn E j e to 
the set of cx-measures does not exceed 16/m, and so, proposition [2] is proved. 

2. The rest of this subsection is dedicated to the detailed description of 
the sequence of hollows (Q £ , I £ ) and to the proof of property (P) for them. 

First consider an auxiliary construction. Take two different points F and 
F' above the line I = {x 2 = 0}, with \OF\ = 2 = \OF% Denote by $ and 
$' the angles the rays OF and OF\ respectively, form with the vector (0, 1). 
The angles are counted clockwise from (0,1). Thus, F = 2(sin $, cos $) 
F' = 2(sin cos$'). Assume, for further convenience, that F is situated 
on the left of F'; thus, one has — tt/2 <$<$'< tt/2. (The case where 
F is situated on the right of F' is completely similar.) Select three positive 
numbers A, A', and 5, and define two ellipses £ and £' and two parabolas 
V and V\ The first ellipse has the foci O and F, the length of its large 
semiaxis is \/l + A, of the small semiaxis, V% and the focal distance equals 
2. The second ellipse has the foci O and F\ the lengths of its large and 
small semiaxes are a/1 + A' and v^A 7 , respectively, and the focal distance is 
also 2. The parabolas V and V have the foci F and F\ respectively, the 
common axis FF\ and the same focal distance 5. Thus, the parabolas are 
symmetric to each other with respect to the bisectrix of the triangle OFF\ 
The parameter S is chosen sufficiently small, so that the point O lies in the 
exterior of both parabolas. See the figure below. 
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In what follows, we shall distinguish between the billiard and pseudo- 
billiard dynamics. The pseudo-billiard dynamics is defined as follows. A 
particle starts at a point (£,0) G I and moves with a velocity (sin (p, cosy?) 
until it reflects from the interior side of S. (Before the reflection it can 
intersect other curves S\ V, V\ or even intersect £ from the outer side, 
without changing the velocity.) Then it moves again with constant velocity, 
until it reflects from the interior side of V. Then, in the same way, it reflects 
from the interior side of V and then from the interior side of £\ and finally, 
intersects I from above to below. Denote by (£', 0) the point of intersection, 
and by — (sin<^', cos the velocity at this point. 

Consider the admissible set: the set of 7-uples (<p, A, A', 5) such 

that all the indicated reflections occur in the prescribed order. This set 
is open and nonempty. Indeed, let 5($, $') be the least of the values 5 
such that one of the parabolas (in fact, both of them simultaneously) passes 
through O. Put cp — $, £ = 0, and take arbitrary values A > 0, A' > 0, 
— 7r/2 <$<$'< 7r/2, < 5 < 5($, $'). The particle with initial data 
= £ = first passes along the large semiaxis of £, then reflects from 
£ , returns along the same semiaxis and reflects from V. Then it moves with 
the velocity parallel to FF\ reflects from V\ moves the large semiaxis of 
the ellipse £\ reflects from it and returns to O along the same semiaxis. 
Thus, the admissible set is nonempty. Under a small perturbation of the 
parameters ip, £, $, A, A', 5, all the reflections are maintained and the 
order of reflections remains the same. This implies that the admissible set is 
open. 

This description determines the mapping <p' = (p J (<p, A, A', 5), 

£' — ^'(v^ £5 ^5 ®\ K from the admissible set to R 2 . This mapping 

is infinitely differentiable. For ip = $ and £ = one has 

^'($,0, A, A', 5) = $'. (11) 

For £ = with arbitrary <p one has 

£'(¥?, 0,$,$', A, A', 5) = 0, (12) 

and 

</?'(</?, 0, A, A', 8) does not depend on 5. 

Indeed, a particle starting at O, after the reflection from £ passes through 
F, then after reflecting from V moves in parallel to FF\ after the reflection 

2 Note that throughout this paper the sign ' (prime) never means derivation. 
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from V passes through F\ and finally, after the reflection from £' returns 
to O (see the figure below). The initial and final velocity of the particle are, 
respectively, (sin ip, cos <p) and — (sin</?', cos <£>'). Denoting by a and a' the 
angles the second and fourth segments of the (5-segment) trajectory form, 
respectively, with OF and OF\ one has a = a\ The angle a is a function of 
p, and </?' is a function of a'; these functions depend only on the parameters 
of the ellipses S and respectively, and do not depend on the parameter 5 
determining the shape of parabolas. 




Using properties of ellipses, one derives the formulas connecting ip, a, 
and <p' = p , (p,0,<Z>,®\\,\ , ,5): 



sm((p— <E>) 



A sin a 



2 + A - 2cosa v / l + A 
It follows that 



, sin(y2'-$') 



A' sin a 



2 + A' - 2cosoVl + A' 
(13) 



dp 1 
dp 



x/A 7 1 + y/X 



VP = $ 



1 + v/A 7 \/A 



(14) 



With fixed $, A, A', and 5 the mapping </?'(</?,£), ^'(v^jO preserves the 
measure, cos pdpdl; = cos p 1 dp 1 d^\ hence 



cos 93 = ±cos<^' 



9y' <9y?' 
dp d£ 
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Using (TT2"j) . one gets that f^-_L_ = 0, hence 



d(p' dip' 

dip S£ 

at d£ 

dip .95 



5=0 



d<p' d? 



5=0 



therefore 



cosy? = ± cos v? t; 

ay? a£ 



(15) 



Putting ip = $, £ = 0, and taking into account (fTT]) . (|T4j) and (fT5l) . one gets 

2 



cos$ = ±cos$' 




(16) 



ip = $ 
5 = 



Define the positive continuous functions A($') and A'($) by the relations 

2 / \ 2 



i + Va 



- cos$', 



l + v 7 ! 7 



- cos$, 
2 



then one has 



ip = <I>; A = A(<J>') 

5 = 0; A' = A'($) 



1. 



(17) 



(18) 



Additionally, taking into account (tHJ) and ( |T7I) . one gets 

cos<3>' <9y?' 



cos $ dip 



■1. 



= $; A = A(<I>') 

5 = 0; A' = A'($) 



(19) 



Recall that y?' = y?'^ 0; ^> ^\ A, A'), that is, the restriction of the function 
y?' to the subspace £ = 0, does not depend on 5. Hence the function ^J^_ Q 

and, by formula (fl5|) . the function f|-L =0 also do not depend on 5. Put 
$ = arcsin(l - 2/m), so that J 1 = [-tt/2, -$ ] and J m = [$„, 7r/2], 
and put A$ = 2/m. The set {($, 0, $, A($'), A'($)) : -$ <$,$'< 
$0) $' — $ > A$} is compact and belongs to the (open) domain of the 
function y?'. Choose a sufficiently large integer value A; = k(e), so that for 



siny? — sin$| < 2 /(km), £ = 0, 
$'_$>A$, A = A($'), 



A' = A'($) 



(20) 
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holds true 

COS$' . n . . . 

1 + £ > 1+£ - 21 

cos$ dip 

Formulas ( 12TI) and (1TTI) mean that under the conditions (l20i) . <y?' is also close 
to Increasing k if necessary, ensure (under the same conditions) that 

cos $' cos 



cos$ cos<^ : 



G [(1 + ^r 1 , l+e]. (22) 



Taking into account ( fT5l) . (|2B . and ( J22l) . one obtains that under the condi- 
tions fl20l) holds true 



9f 
<9£ 



G[(l + .)- 2 , (l + e) 2 ]. (23) 



3. Now we proceed to the description of the hollow (fl e , I e ). 

(a) If 2 < i 7^ a(i) < m — 1, divide the interval J* into /c subintervals 
J 4 '- 7 of equal measure A, going in increasing order: J % = Uj =1 J l '\ A(J M ) = 
2/ (km) for any j = l,...,k. Recall that d\(ip) = cosipdip and the value 
k = k(e) is defined above. Without loss of generality assume that k(e) — > oo 
as £ — > 0. 

To each pair of intervals, J w and J t7 W ,J °, we apply the construction de- 
scribed above, see fig. 9. Namely, draw arcs of ellipses =AB, £\j =A'B' 
and arcs of parabolas Vij, Without loss of generality suppose that 

i < o~(i). The angles AOB and A OB' correspond to the angular inter- 
vals J M and J^w^', respectively. The foci F = Fij and F' = F\ t j belong 
to the intervals OA and OA\ respectively. The endpoints of the arcs Vi.j 
and V\j also belong to the intervals OA and OA\ respectively. The angle 
corresponding to the ray OA (and therefore to the left endpoint of the in- 
terval J* J ) will be denoted by $ = and the angle corresponding to the 
ray OA' (and therefore to the right endpoint of the interval J "^') will be 
denoted by $' = Denote A = Ay := A($') and A' = A'y := A'($), ac- 
cording to the formula (j!7p . Next, select a value 5 = Sij and draw two 
curves (lateral reflectors) in such a way that (i) each of the curves con- 
tains an arc of parabola (the first curve contains Vij, and the second one, 
V\j), an arc of circumference centered at O, and three radial segments; (ii) 
these curves do not intersect the intervals whose endpoints belong to the set 
{F a j3, Fly t s : (a,j3) ^ (i,j), (j,b~) ^ j)}' this will guarantee free pas- 
sage of particles from one parabola to another; and (iii) the A-measure of the 



24 



angular interval occupied by each lateral reflector does not exceed e/(km). 
On the figure below, the angular reflectors are the curves joining the points 
A and C, and the points A 1 and C . 




Notice that -$ $o and $' - $ > A$. Indeed, $ and $' do 

not belong to the intervals J 1 = [— 7r/2, -$ ] and ^ m = [$o> tt/2]. On the 
other hand, one has $'-<!> sin<l' - sin$ = A([<§, $']) > 2/m = A$. 

Introduce the shorthand notation (p'((p, £) = £, Ajj, AYj, d>ij). 

According to fl2TT) and fl23|) . for (/? G J*' J holds true 



cos$' . r/1 , i 

^7r^'°) G t( 1 + £ ) ' 1 + 

cos <P ay? 



and 



<9£' 



g [(l + ^r 2 , (i+ e ; 



(24) 



(25) 



According to ffTTT) . one has y?'^ 0) = ^'i this equality and the formula ([2 
imply that for ip G and y?' = y?'(y2, 0) on e has 



cos$' tp' - $' 



cos $ y? — $ 
On the other hand, one has 



G 1 + e]. 



(26) 



cos$|/' J | = — (1 + o(l)), 
/cm 



(27) 
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cos$'|J CT(i) ' j | = — (l + o(l)), (28) 
Km 

with o(l) being uniformly small over all i, j as e —>■ 0, and |J| being the 
Lebesgue measure of J. (Recall that the parameters $, and the inter- 

vals J M implicitly depend on e.) 

Choose closed intervals J M C and J a< ^j c J^W^ satisfying the 
following conditions: (i) <^'(J J,jf x {0}) = J <7 W>- J ; (ii) some neighborhoods of 
J JJ and j^W'- 7 belong to J M and J ™'- 7 , respectively; and (iii) the pseudo- 
billiard trajectory with the initial data (</?, 0), <p G J M does not intersect 
the neighbor lateral reflectors (that is, the lateral reflectors corresponding to 
the intervals J^ +1 and J"®*- 1 , if j ^ 1, k; if j = 1 then J'Wj-i should 
be replaced with j^W -1 ^ and if j = k then J M+1 should be replaced with 
Note in this regard that the neighbor lateral reflectors occupy a small 
part of the angular intervals J M and J a ^>i (represented on the figure by the 
arcs AB and B'A'). Other lateral reflectors will not be intersected, by the 
choice of lateral reflectors. 

By virtue of (1261) . (1271) . (1281 and because of smallness of the angular 
intervals occupied by the lateral reflectors, J J,jf and J a ^j may be chosen in 
such a way that the ratios A( J iJ )/X( J iJ ) and A( J a{i) ' j )/\{ J a ^) uniformly 
(with respect to i, j) tend to 1 as e — > 0. Thus, a billiard particle going 
from O in a direction ip G J , J , makes the same reflections and in the same 
order as under the pseudo-billiard dynamics: first, reflection from £ij, then 
from Vij, from V\j, from S\j, and finally, the particle goes back to O in 
the direction y?'^ ) e J ctWj '. 

Choose Ojj in such a way that the following conditions are fulfilled: if 
(£, v 2 ) ^ [— Qjj] x then (i) the corresponding billiard trajectory does 
not intersect the lateral reflectors and the indicated order of reflections is 
preserved; (ii) p'foO G J^ j ; (iii) §(<M)| e [(1 + ^)- 3 , (1 + e) 3 }. 
Analogously, choose a a (i),j m such a way that the conditions are fulfilled: 
if (£, v?) G [— a<r(i),j, a-a(i)j] x J "'- 4 '''- 7 then (i) he billiard trajectory does not 
intersect the lateral reflectors and the order of its reflections is reverse; (ii) 

p'(^5)eJ«; (iii) f|(<M) e[(l + e) -3 , (1 + e) 3 }. Note that the values 
dij and a<T<i),j implicitly depend on e. 

Select a e < min^- in such a way that a e — > as e — > and denote 
I £ = {-a £ , a £ ) x {0}, I e = (-a £ (l + £)~ 3 , a £ (l + ep 3 ) x {0}, and J\ = 
J 1 := Uj J*'- 7 . The part of the boundary of f2 e related to the angular intervals 
J*'- 7 and J a w>^ under consideration is formed by the arcs of ellipses Sij, 
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£\j and the corresponding lateral reflectors. Then a billiard particle with 
initial conditions (£, </?) G 7 e X J 4 '- 7 after making four reflections will intersect 
Z at a point (£', 0) G J E , and the angle at the point of intersection will be 
<Pa,,i.(£> ( P) = e J CT(i),j C J CT(i) . Thus, the set of values (f, G J £ x J* 
such that <PQ StIe (£,,<p) & J a ^ is contained in the set (I £ x J' 1 ) \ (j £ x J*J, 
whose measure is vanishing as e — * 0. Q.E.D. 

(b) If 2 < z = er(«) < m — 1, the corresponding part of the boundary 
is the arc of circumference of radius 2 with the center at O occupying the 
angular interval J\ that is, the set {2(siny9, cosy?), <p G J 1 }. Next we will 
show that for all values (£, 0) G I E , <p G J*, except for a portion of order 
o(l), the corresponding billiard particle makes one reflection from the arc 
and then goes back to I £ in the direction </?' G J J . 

For all values ^ 6 J 1 , except for the union of two intervals of vanishing 
length (each of the intervals is contained in J 1 , has the length 2 arctan(a E /4), 
and contains an endpoint of J 1 ), the particle starting at (£, 0) G I £ in the 
direction (p will reflect from the indicated arc of circumference. Let ip G J 1 be 
the angular coordinate of the reflection point. By (£', 0) denote the point at 
which the reflected particle intersects the straight line I. One easily verifies 
that 

7 + 77 = cos ^- (29) 
Ifl < a e , (30) 



One has 
hence 



1 I / 1 



> (31) 



From ([29"]) it follows that |£ + £'|/l££'l = I cos^l < 1, and taking into account 
( 130]) and (131]) . one finds that |£ + £'| < a 2 £ / {1 — a £ ). This implies that for 
all values (£, 0) G I £ , except for a set of measure 0(a £ ), the second point 
of intersection of the billiard trajectory belongs to I £ , moreover the velocity 
at this point, ^ 7 (£, <p), belongs to A/2arctan(o s /4)(^ 1 )' the neighborhood of 
J 1 of radius 2 arctan(a e /4). This finally implies that for all (£, <p) G I £ x J*, 
except for a portion of order 0(a £ ), holds <Pq j (£, y?) G J*. 

(c) The parts of the hollow's boundary, corresponding to J 1 and J m , 
are formed by smooth curves joining the corresponding endpoints of I £ and 
the points 2(sin <3> , — cos $ ) an d 2(sin $ , cos <& ), respectively. The unique 
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condition on these curves is that they can be parametrized by the monotoni- 
cally increasing angular coordinate. For those values cr(l), a(m) that do not 
coincide with neither 1 nor m take just the arcs of circumference of radius 2 
corresponding to the angular intervals J a( - 1 \ J°"( m ). 

Consider the union of all the elliptic arcs Sij, Slj introduced in item (a), 
all the arcs of circumference defined in items (a) and (b), and the two curves 
introduced in this item (c). Let us call this union the main element. Each 
lateral reflector is a curve; select it in such a way that both its endpoints 
belong to the main element. Finally, the curve dfl £ \ I £ is the union of all the 
lateral reflectors and the part of the main element visible from O (that is, 
which is not shielded by the adjacent lateral reflectors). Thus, the definition 
of the hollow (Q e , I e ) is complete. 

On the figure below, there is shown a particular hollow (Q £ ,I £ ) corre- 
sponding to the permutation a — ( \ \ \ \ \ ) . The angular intervals 

J 1 , . . . , J 5 are separated by dotted lines. The family of hollows (Q e , I e ), with 
vanishingly small e, has the following property: for almost all particles with 
the initial direction from J 2 (resp. J 3 , J 4 ), the final direction will belong to 
J 4 (resp. J 3 , J 2 ). On the figure, there is shown the trajectory of a particle 
with the initial direction tp e J 2 and the final direction (p + e J 4 . The par- 
ticle makes a reflection from an elliptic arc, then two reflections from (very 
small) parabolic arcs, and finally, again from an elliptic arc. According to 
our notation, these arcs are £2,2, ^2,2, ^2,2, and £\,2- 
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4 Concluding remarks and applications 

Physical bodies in the real world have atomic structure and therefore are 
disconnected. This is a reason for using (generally) disconnected sets Q m in 
the definition of a rough body. In future we intend to turn to propose and 
study the notion of a three-dimensional rough body, where the connectivity 
assumption is absolutely useless; this is another reason. By removing this 
assumption, the consideration in two dimensions (namely, proof of lemma 
[T]) is made somewhat more difficult, but at the same time prerequisites for 
passing to the three-dimensional case are created. 

In fact, the notions of "disconnected" (as everywhere in this paper) and 
"connected" rough bodies are equivalent. There is a natural one-to-one corre- 
spondence between the equivalence classes in the connected and disconnected 
casesjl the former classes being subclasses of the latter ones under this cor- 
respondence. 

Let us now consider applications of theorem to problems of the body 
of minimal or maximal aerodynamic resistance. A two-dimensional convex 
body B moves, at constant velocity, through a rarefied homogeneous medium 
in M 2 , and at the same time slowly rotates. The rotation is generally non- 

3 more precisely, equivalence classes formed by sequences of connected / disconnected 
sets 
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uniform; we assume that during a sufficiently long observation period, in a 
reference system connected with the body the body's velocity is distributed 
in S* 1 according to a given density function p, with J gl p(v) dv = 1. The 
medium particles do not mutually interact, and collisions of the particles 
with the body are absolutely elastic. The resistance of the medium to the 
motion of the body is a vector-valued function of time. After averaging it 
over a sufficiently long period of time, one gets a vector. We are interested 
in the projection of this vector on the direction of motion; for the sake of 
brevity, it will be called mean resistance, or just resistance. The problem is: 
given B, determine the roughness on it in such a way that main resistance 
of the resulting rough body is minimal or maximal. 

A prototype of such a mechanical system is an artificial satellite of the 
Earth on relatively low altitudes (100 -j- 200 km), with restricted capacity of 
rotation angle control. The satellite's motion is slowing down by the rest of 
atmosphere; the problem is minimize or maximize the effect of slowing down. 
The problems of resistance maximization may also arise when considering 
solar sail: a spacecraft driven by the pressure of solar photons. 

The initial velocity of an incident particle (in the reference system con- 
nected with the body) is —v, and the final velocity is t> + ; therefore, the 
momentum transmitted by the particle to the body is v + v + . The projec- 
tion of the transmitted momentum on the direction of motion of the body 
equals 1 + (v, v + ). Averaging this value over all particles incident on the 
body within a sufficiently long time interval, one gets the mean resistance. 
The averaging amounts to integration over p{v) du B (v, v + , n); that is, mean 
resistance of the rough body equals 

RM = JJJ^(l + {v,v + ))p(v)du B (v,v + ,n). 

Using theorem 1 and Fubini's theorem, one rewrites this formula in the form 

R(u B ) = I dr B (n) II (l + {v,v + ))p(v)dr] Btn (v,v + ), (32) 
Js 1 J Jr 2 

where r] BjJl e A n . Thus, the minimization problem for R{v B ) reduces to 
minimization, for any n, of the functional JJ T2 (1 + (v, v + )) p(v ) dr](v , v + ) 
over all r\ e A n . Using the notation introduced in subsection 3.3, one comes 
to the problem: 

inf // (l + cos((p-(p + ))g(ip)dri(ip,(p + ), (33) 
^ A JJn 



30 



where g(p) = p(v) for p = Argt> — Argn. This problem, in turn, by sym- 
metrization of the cost function reduces to a particular Monge-Kantorovich 
problem: 



where c(ip, ip + ) = (1 + cos(<p — <p + )) e 2 g an d A A A is the set of measures 
77 on □ having both marginal measures equal to A: tt^t] = A = ir^+f]- Recall 
that A is defined by d\(v) = cos if dip. 

The problem (1341) can be exactly solved in several particular cases. Con- 
sider the case of uniform motion, where the function p, and therefore g, is 
constant, and thus, one can take c(ip, ip + ) — | (1 + cos(<p — y? + ))0 Note that 
jF(r/ ) = 1 and therefore resistance of the smooth body is equal to its perime- 
ter: R{vb) = fgi c?Te(n) ^(rjo) = \dB\. (Recall that the measure rj belongs 
to A and is supported on the diagonal ip + = —ip.) The minimization problem 
( 134)) for constant g was solved in [H]: one has infe_R(z/e) = 0.9878... • \dB\, 
the infimum being taken over all roughenings of B. 

Note that the corresponding maximization problem for (IMI) has the trivial 
solution, which does not depend on the function g: 77 = 77*, the measure 77* G 
A being supported on the diagonal ip + = ip. One has sup B R(u B ) = k \ dB\, 

where k = yJl^Qi^) cos pdp^j j (^fl^^if) cos 3 pdp^ > 1; in the case 
of uniform rotation one has k = 1.5. The maximization problem was studied 
in more detail in [20] . 



The construction is simple (see the figure), but its description is a bit cum- 
bersome. 

Take a point in the interior of B and connect it by segments with all 
vertices. The polygon is thus divided into several triangles; fix i and m and 
consider the triangle with the base b,, the ith side of B. Denote by d^Q™ 1 ) the 
diameter of the orthogonal projection of f2™ on the straight line containing 
/™; one obviously has rf(f2™") > Fix a positive number k < |/ J m |/<i(f2™'). 
Take a rectangle II 1 contained in the triangle and such that one side of II 1 

The normalization constant 3/8 is taken for further convenience. 



inf T{rf) , where 




(34) 



Appendix A 
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belongs to b,. By 81 denote the total length of the part of bj which is not 
occupied by this side. 

For the sake of brevity, the image of a set under the composition of a 
homothety with positive ratio and a translation will be called a copy of this 
set. Take several copies of f2™ (copies of first order) that do not mutually 
interact, belong to II 1 , the corresponding copies of If 1 belong to bj, and the 
portion of the side of II 1 occupied by them is more than k. 



Next, take several rectangles that do not mutually intersect and do not 
intersect with the chosen copies of belong to IT 1 , and have one side 
contained in bj. Denote by II 2 the union of these rectangles and by 82, 
the total length of the part of the side of II 1 which is not occupied by the 
rectangles from II 2 and by the copies of J™. Next, for each rectangle from IT 2 
choose several copies of Q™ (copies of second order) in the way completely 
similar to the described above (see fig. ??). 

Continuing this process, one obtains a sequence II 1 , II 2 , ... of unions of 
rectangles and collections of copies of Q" 1 of 1st, 2nd, . . . order. Choose the 
rectangles in such a way that S 1 + 5 2 + ■ ■ ■ < l/m and Area(II 1 ) < l/m. 
Finally, choose k such that the total length of sides of rectangles from Tl k+1 
contained in b« is less than l/m, and take the collection of copies of Q™ 1 of 
order 1, 2, . . . , k (we shall call it full collection). The total length of the part 
of bi not occupied by the corresponding copies of i|™ is less than 2/m, and 
therefore goes to zero as m — > 00. 



hi 
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By definition, the desired set Q m is B minus the union of full collections 
of copies of Q™ over all i. 



Appendix B 

We prove here slightly more than needed. 

Statement 3. Let A = (a>ij)ij=i be a symmetric matrix, with being 
nonnegative integers. Denote ni = Yl*j=i a ij- Then there exist matrices 
Bij — {bij)n,v of size rii x rij such that e {0, 1}, Bjj = Bj i} the sum 
of elements in Bij equals a^, and the block matrix D = (Bij) contains ex- 
actly one unit in each row and each column. 

Note that for some values i — i±, i 2 , ■ ■ ■ it may happen that rij = 0, that 
is, — for all j — 1, . . . , k. Then the corresponding matrices B^ have the 
size x rij, that is, are empty. In this case D coincides with the block matrix 
D' = (B^) having the rows i±, ■ ■ ■ and columns i±, i 2 , ■ ■ ■ crossed out. 



Proof. The proof is by induction on k. Let the statement be true for k — 1; 

l ij)ij=2'i 



prove it for k. Take the matrix A = (ajj)*,- =2 ; there exists a block matrix 



B = (Bij)ij =2 satisfying the statement. Note that the order of B^ is rij x -hj, 

where hi = Y^=2 a ij = n i ~ a n- Define the matrices By as follows. 

(a) Put B u = diag 0, . . . , 0}. 

an 

(b) Put ^ ail+1,1 = = ^ a n+ a i2i a l2 _ ^. ^ail+dl2+l,l _ _ ^an+ai2+ai3 5 ai3 

1; . . . ; 6 «ii+-+ fl i.fc-i+ 1 ' 1 = . . . = 6 «ii+-+«ifc.«i* = i ; the other elements of the 
matrices B^, j = 2, . . . , k are zeros. Thus, on the diagonal of B\j starting 
from the element at the first column and the (an + a 12 + . . .+ai J _i + l)th row, 
the first a\j elements equal 1, and the remaining elements on this diagonal 
and all the elements off the diagonal are zeros. This defines the matrices Bij, 
j = 2, . . . , k. The matrices B a , i = 2, . . . , k are determined by the condition 
Bn = B^. 

(c) For % > 2, j > 2 define the matrix B^ as follows. For fi < an or 
v < aij, put bij = 0, and for /j > an + 1, v> aij + 1, put = au ' u aiJ . 
Thus, in the obtained matrix B^, the right lower corner coincides with the 
matrix B^, and all the remaining elements are equal to zero. The number of 
rows of this matrix equals an + hi = n iy and the number of columns equals 
a^ + hj = nj. One obviously has Bj- = Bji. 
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One easily verifies that bfj = and that each row and each column 
of the obtained block matrix D = (-B^)* - =1 contains precisely one unit. 

□ 
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